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Abstract 



The dynamical structure of an atomic Bose-Einstein condensate limits 
the efficiency of the condensate in cooling slow impurity atoms. To illustrate 
the point, we show that an impurity atom moving in a homogeneous zero- 
temperature condensate is not scattered incoherently if its velocity is lower 
than the condensate sound velocity c, limiting cooling to velocities v > c. This 
striking effect is an expression of superfluidity and provides a direct means 
to detect the fundamental property of superfluidity in atomic condensates. 
Furthermore, we show that the fermionic lithium-isotope, ^Li, is a reasonable 

candidate for sympathetic cooling by a ^'^Na-condensate. 
PACS numbers(s):03.75.Fi, 05.30.Jp, 32.80Pj, 67.90.-hz 
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The significance of tlie observation of Bose-Einstein condensation in alkali-atom traps 
||T|] can hardly be overstated : these experiments represent the first explorations into the 
realm of dilute atomic quantum gases. It is thus unfortunate that the same technology is 
not expected to cool fermionic atoms into the degenerate regime. The problem stems from 
the antisymmetry which forbids s-wave scattering of indistinguishable fermions, implying a 
vanishing thermalization rate in evaporative cooling at low temperatures 0. To overcome 
this difficulty, the JILA-coUaboration proposed an interesting alternative P] : the fermions 
can be cooled by bringing them in contact with an atomic condensate. This scheme of 
sympathetic cooling would be an interesting and important application of BEC, cooling not 
only fermionic alkali atoms, but any gas that cannot be cooled by means of the traditional 
cooling techniques, such as inertial gas atoms or perhaps even molecules 0] . 

In this letter, we discuss the efficiency of such cooling schemes 0. Whereas one generally 
expects sympathetic cooling to be most efficient with a cooling system of lowest possible 
temperature, we point out that this is not necessarily true for a Bose condensate. In fact, we 
show that for a zero-temperature condensate, the cooling of a low-density gas of 'impurity 
atoms' (distinguishable from the condensate atoms) exhibits a limit due to the superfluidity 
of the condensate: impurity atoms slower than the condensate sound velocity are not scat- 
tered incoherently and cannot be cooled down to lower velocities. At finite temperatures the 
condensate contains a normal fluid fraction and the impurity atoms cool down further, but 
at a reduced rate. This aspect of sympathetic cooling provides the opportunity to observe 
the fundamental property of superfluidity. 

At zero temperature the condensate scatters impurity atoms by creating quasi-particles. 
In a condensate of bosons with mass and density hb, the energy dispersion of the 
quasi-particles is Ug = cq^l + (q/kc)'^, where = v^lGvm^a is the inverse coherence length 
and c = kc/^rriB the condensate sound velocity. In the long wavelength limit, q <C k^, 
the momentum region of importance in describing slow impurity atom scattering, ujq ~ cq 
is phonon-like. We can then repeat Landau's argument for the absence of resistance in 
superfluid flow, which translates here into the absence of incoherent scattering for impurity 
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atoms moving at a velocity v < c. The argument is based on the observation that energy 



energy transfer AE = cJ2n l^tnl, where the q„- vectors represent the momenta of the phonons 
created in the scattering, giving a total momentum transfer q = J2n^.n- From the energy 
difference experienced by the scattered impurity atom of mass rrij, AE = v-q— g^/2m/ < vq 
and since q < J2n lQn| = AE/c, we find that AE < {v/c)AE, which cannot be satisfied if 
V < c. 

To obtain a quantitative understanding of BEC sympathetic cooling, we discuss a some- 
what idealized model of a homogeneous zero-temperature BEC, cooling an impurity-atom 
system of such low density that 1. the interactions between impurity atoms can be neglected, 
2. the impurity system remains non-degenerate throughout the cooling process and 3. we 
can neglect the heating of the condensate. The occupation number nk of the plane wave 
state of momentum k evolves in time according to a simple Boltzmann equation : 



where w(kj„; kj) is the rate for scattering of an impurity atom from an initial momentum kj„ 
to momentum kj. The w-rate is proportional to the dynamical structure of the scattering 
(or cooling) system. We derive this result for an impurity atom that interacts weakly with 
the bosons by means of a pseudo-potential, H' = \ Y^f=i5{r — rj), where r represents the 
impurity atom position and the position of the i-th boson, and where A = 27ra/B/yU, with 
ajB the scattering length for the impurity-boson scattering and fi the reduced mass, = 
m]^^+mj^. The interaction matrix element between the initial state, |\l/o,kj„) = |\l'o)(8'|kj„), 
where |\l/o) denotes the BEC ground-state, and the final state kj) = ® |kj), is then 
equal to 



and momentum cannot be simultaneously conserved if f < c. To see this 0, consider the 
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where q is the momentum transfer km — k/, pq the Fourier transform of the particle density 
and V the macroscopic volume of the homogeneous system. Representing the initial and 
final state energies by Eq + e(kj„) and Ef + e(k/), we obtain for the scattering rate in a 
Fermi-Golden rule calculation: 

«;(k,„;k^) = 27r^ |(^^,k;|if'|^o,k..)|' 

I*/) 

X 5(E^ + e(k/)-[Eo + e(k,„)]) 
= ^^{y) 5(q, [e(k„) - e(k;)]) , (3) 

where we introduced the dynamical structure factor S'(q,co') = J2\^s>f) l(^/|Pql^o)P^(-E'/ — 
Eq — uj). In the thermodynamic limit, X^q ^ ^/(2vr)^/ d^q, the Boltzmann-equation (|I]) 
takes on the following form : 

2 

I - / ■s(q, V ■ q - q'^/2mi)nk{t) 
+ j (fq s(q, V ■ q + q^ /2mi) nk+q(t) , (4) 



where s denotes the dynamical structure factor density, s(q, cj) = [S{c[,uj) jV]. For a zero- 
temperature dilute condensate s is approximated accurately by 

s(q, u) = ub 5{uj - uj^) , (5) 

describing the creation of a single quasi-particle of momentum q and energy ujq (defined 
above). Although the above derivation assumes weak interaction between the impurity and 
the condensate (reasonable if w > c), the rigorous derivation of the quantum Boltzmann 
equation ^ shows that Eq.(||) has a broader range of validity. 

With Eqs.(^ and (|^), we calculate the rate dQ/dt = — J d^ke(k)hi^ at which energy is 
transferred from the impurity atoms to the condensate. 



^ = j^ld'kn^ /rf'g^.^(q,Kk)-e(k-q)]), (6) 



where is normalized so that / d^k nk equals the number of impurity atoms. Performing 
the remaining q-integration, we find 
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where Fcik) = if f < c, and 
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if f > c. A classical system of noninteracting particles of mass has a dynamical structure 
factor density s(q, cj) = ns^ioJ — jlrtiB)-, giving the above result with Fc 1. Thus the 
Fc{v) -function, shown in Fig.l, corrects for the condensate dynamical structure. At f ^ c, 
F(.[v) ~ 1 and the energy transfer rate for cooling fast impurity atoms by a condensate is 
the same as the rate for cooling by noninteracting atoms at rest. At f > c, the condensate 
dynamical structure factor decreases the rate of energy transfer until at f < c the impurity 
atoms are not scattered and no transfer of energy takes place, confirming the picture derived 
from the Landau argument above. 

Even if interactions between the impurities cannot be neglected (giving impurity- 
impurity collision terms in the Boltzmann equation), the heat transfer rate of Eq.(|^) is 
correct, provided the impurity system is non-degenerate and interacts weakly with the con- 
densate. In the limit of strong impurity-impurity interactions, these interactions 'ther- 
malize' the distribution during the cooling process and n\Sf) remains Maxwellian, 
n^Sfy ~ exp \—k'^l2mikBT(i)\^ where ks denotes the Boltzmann constant. We can then cal- 
culate the rate of temperature decrease in accordance with the heat transfer rate, dT/dt ~ 
— [dQ/dt]/C, where C is the heat capacity, C = j d^k nke(k) = / d^k t'^i\^)n\^/kBT'^ ■ 
This gives dT/dt ^ —T/tc^ where the cooling rate depends on the thermal velocity, 
Vt = \/2mBkBT, as t^^ = {ii^ /miniB) (IG/lS-y/vr) x AnnBoj^VT Fs{vt) ~ {fi'^ /mjuiB) x 
AimBO^jBVT Fs{vt), where the Fs-factor corrects for the condensate structure factor, express- 
ing the corresponding suppression of the cooling rate, Fs = J d^k n]Jx,^Fc{k) / J d^k rikfc^ 
(see Fig.l). 

However, the interaction with a dense condensate generally modifies the condensate. 



We return, therefore, to to the low-density impurity system where nk(t) does not remain 
Maxwelhan (although the above cooling rate is still a reasonable estimate) and we solve the 
Boltzmann equation numerically. In Fig.2, we show the time evolution of an isotropic im- 
purity distribution that is cooled by a zero-temperature condensate. The distribution starts 
out as a Maxwellian with vt = 3c, and progressively narrows in velocity-space but does not 
converge to zero velocity - instead the position of its maximum converges to a value slightly 
below c: the signature of superfluidity. The distributions are plotted at times t = OAnr, 
n = 0, ■ ■ ■ , 10, where is the cooling rate for vt = c, = Ana'j^nBcl^'^ /mimB). The 
mass-ratio of the impurity and condensate atoms in the calculation is 6/23, corresponding 
to the fermionic ^Li-isotope cooled by a ^^Na-condensate. 

The choice of the elements and their atomic states is crucial in sympathetic cooling since 
processes such as spin-flip scattering limit the lifetime of the impurity system in contact 
with the condensate. ^^Na, the atom used in the high-density BEC-systems created at MIT, 
and ^Li are promising choices: their atomic properties are well understood and the spin-flip 
rates relatively small. 

The magnetic traps in which EEC has been realized confine only a fraction of the hy- 
perfine states, namely the low magnetic field seekers. At small field intensities, the low field 
seeker states are (/ m) = (| |), (| |) and (| -i) for "^Li, and (2 2), (2 1) and (1 -1) for 
^^Na. In the elastic approximation [0, also known as the Degenerate Internal States (DIS) 
approximation |Tl|], one describes the scattering of atoms in these hyperfine states in terms 



of the singlet and triplet scattering amplitudes. At the relevant low energies, the scattering 
is dominated by the s-wave contribution and the triplet (singlet) phase shift St{s) can be 
represented by —kaT(s)i where k is the relative momentum and aT(5) is the scattering length 
for the triplet (singlet) scattering (see Table |). The scattering length for elastic processes 
(responsible for cooling and thermalisation) is then a ^ Pso-s + Pto-t where Ps and Pt 
represent the probability that the scattering product is in the singlet and triplet state. For 
example, ^^Na-^^Na scattering in (1 -!) + (! -1) gives a ~ + jq^s = 69ao, in good 



agreement with SGtga'^o HH and close to 52 ± Sag |jT3|. The spin- flip cross section is 



(Js/ = Afs/ vr(aT - as)^ (9) 

where M^f is a coefficient that depends on the specific states involved in the process. This 
multiphcative factor vanishes if the initial state of the scatterers is of a pure triplet or singlet 
character, implying that it is possible to choose the atomic states so as to avoid spin-fiip 
processes all together. In cases where spin-fiip occurs, their rate may still be reduced if the 
singlet and triplet scattering lengths are nearly equal, as for the ^^Rb-condensates observed 
at JILA [i|l9|,|0|. 

For ^Li -®Li collisions, only collisions between (/ m) states (| |) + (| |), (| |) + (| 
and (| -|) + (| -|) are without spin-fiip. Because of the fermionic nature of ^Li, these are 
all forbidden to s-wave collisions (no cooling possible). For ^^Na-^^Na collisions, the only 
scattering products not leading to spin-fiip decay are (2 2) + (2 2) and (1 -1) + (1 -1). The 
collision channel (2 2) + (2 2) is purely elastic with a given by aj- in Table |. (1 -1) + (1 -1) is a 
mixture of both singlet and triplet and many exit channels could lead to trap loss. However, 
at low temperature, all these channels are closed. 

For ^'^Na-^Li collisions, the combinations of states without spin- fiip are (| |) + (2 2) and 
(i -|) + (1 -1). Only elastic collisions are possible for (| |) + (2 2) but for -|) + (1 -1), 
many decay channels exist, although not physically accessible at low temperatures. 

The scattering lengths for ^Li reported in Table |I| were derived by Abraham et al. |n 



and those for ^^Na were calculated by Cote and Dalgarno JTSf. They are in agreement with 
recent estimates For the mixed case ^^Na-^Li, we computed the scattering lengths 

using singlet and triplet model potentials formed by the extension of ab initio data with long- 



range tails. The ab initio data originate from Schmidt-Mink et al. ||T6[ and extend to large 
distances (30ao)- Beyond the ab initio region, the potential is described by the dispersion 
coefficients of Marinescu et al. [|17| and takes the form — Ce/r^ — Cg/r^ — Cio/t^° =1= Ar"e~^^. 
The =F sign of the exchange term stands for the singlet and triplet potential respectively. 
The parameters were computed using expressions given by Smirnov and Chibisov |T8| and 



take on the following values (in atomic units): A = 0.0124, a = 4.626 and (3 = 1.2445. The 
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two regions were smoothly joined using a cubic spline fit. 

The spin-flip rate coefficients (if spin-flip scattering occurs and M^/ ~ 1) can be esti- 
mated with the values in Table | by 7^ = (vase) where v is the velocity of the pair of atoms 
undergoing a spin-flip, and where the average is taken over a Maxwellian velocity distribu- 
tion. For ^Li-^Li, TZ is quite large due to its large scattering length, and still sizeable for 
23Na-23Na, but one order of magnitude smaller for ^^Na-^Li (see Table |T|). This situation 
is similar to the case of '^^Rb, where TZ ~ 2.2 x 10~^'^cm^/s and sympathetic cooling of two 
different hyperfine states has been observed 0], presumably because the singlet and triplet 
scattering lengths are nearly equal. This implies that the four other possible collision chan- 
nels for 23Na-6Li, (| i) + (2 2), (i -i) + (2 2), (| §) + (! -1) and (| i) + (l -1), could lead to 
sympathetic cooling, as long as the ^Li states are not mixed (otherwise spin-flip among ^Li 
states will take place). In that case, sympathetic cooling is feasible if the cooling rate 
exceeds the rate of spin-flip scattering ubTZ or 

M.;(^Ia^)7^] « 1, (10) 

\ aiB J ymjuiBj 

where the jb subsripts denote that the scattering lengths are calculated for impurity-boson 
scattering. 

Finally, we remark that the best option for sympathetic cooling of ^Li by ^^Na in the 
(1 -1) state (such as in the MIT condensate) is to trap ®Li in the state (| -|); this ensures 
a small trap loss from spin-flip. On the other hand, a condensate of ^^Na-atoms in the (2 
2)-state with ^Li in the (| |)-state is the best option since no spin-flip is possible. In both 
cases, there is a sizeable scattering length, ajB, hence efficient sympathetic cooling. 

In conclusion, we have shown that the superfluidity of atomic Bose condensates limits 
their efficiency for sympathetic cooling. This effect provides a possibility for direct obser- 
vation of the superfluidity of the dilute atomic trap condensate systems. Furthermore, we 
have discussed the feasibility of such sympathetic cooling, and its dependence on the atomic 
scattering parameters. For the special case of cooling ^Li and ^^Na, we showed that sympa- 
thetic cooling is possible and we discussed the best options for the choice of hyperfine states 
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to be trapped. 
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TABLES 



6Li 23Na 23Na.6Li 



ax -2160 77.3 31.1 

as 45.5 34.9 39.2 



TABLE I. Singlet and triplet scattering lengths in Bohr radii ao, for isotopically pure and 
mixed alkali gases. 

r(K) ^Li-^Li 23]^g^_23]sra 23]^g^_6Li 

lO-'^ 1.1 X 10-9 2.1 X 10-^^ 1.2 X 10-1^ 

10"^ 3.6 X 10-9 6.8 X 10-^3 3.8 x IQ-^'' 

10-^ 1.1 X 10-^ 2.1 X 10-^2 1.2 X 10-^3 



TABLE II. Spin-flip rate coefficients TZ for collisions between ®Li and ^^Na atoms in units of 
cm^/s. 
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Figure Captions 

Fig.l : Plot of the Fc{v) (full line) and Fs{vt) (dotted line) functions that quantify the 
suppression of the BEC cooling efficiency as compared to cooling by atoms at rest (for 
mi/niB — 6/23). Fc{v) represents the suppression of heat transfer and Fs{vt) gives the 
cooling rate suppression of a Maxwellian impurity distribution with thermal velocity vt- 
Fig. 2 : Velocity distribution of the impurity atoms as they are cooled by a zero-temperature 
BEC, plotted at i = OAnT,n = 0,1,..., 10, r = [ATrajsnBcdi^/mimB)]'^- The initial 
distribution is a Maxwellian with vt = 3c. Note that the distribution does not converge to a 
zero-temperature Maxwellian, but the position of its maximum converges to a value slightly 
below c: the signature of superfluidity. 
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